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Methods for determining the dependence of the average distance between particles in monodisperse compos¬ 
ites are analyzed. A more accurate method of determining the average distance between particles in disperse 
systems with an indication of the coordination number taking account of the volume density, connectedness, 
and degree of compaction of the system on the basis of the corpuscular and percolation models is proposed. A 
comparison of experiments and the calculations performed with the more accurate method show good agree¬ 
ment. The method is applicable for composite, porous, and colloidal systems. 
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The structure of composite and porous materials largely 
determines their physical characteristics [1 -5]. The basic 
structural parameters of these systems are particle size and 
interparticle contacts, character and density distribution of 
the solid phase in the disperse medium (matrix). 

The properties of the macro- and microscopically 
inhomogeneous powder materials are described with differ¬ 
ent fractal [1,2] and structural-percolation [3 - 5] models. 

Apparently, fractal models are desirable only for nano¬ 
materials, composites with three (or more) fractions, and 
self-organizing (in which the elements of the universe’s cre¬ 
ator are clearly evident) systems. Structural and percolation 
models are physically more rigorous for porous and compo¬ 
site materials consisting of particles larger than 1 pm. For 
example, it is convenient to use the structural model pre¬ 
sented in [4, 5] for permeable and composite monodisperse 
systems, where indices such as the following are used: D — 
the average particle size in the system (material), P — poro¬ 
sity of the system, V — the relative volume packing density 
of the particles in the system, d — size of the interparticle 
contacts. The missing parameters are easily determined from 
the engineering expressions: Z— average coordination num¬ 
ber; S 0 — specific surface area; r — average radius; E, — the 
sinuosity of the capillaries (for porous materials); N c — num¬ 
ber of contacts; and, S oc — area of the interparticle contacts 
per unit volume. Here we start from the fact that unit volume 
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of the material being modeled contains N= 6(1 - P)/(nD 3 ) 
particles: 

V= 1 -P = 

(1 -2/Z) 2 /[8(1.077- 1/Z-Z- 116 )(0.5- h/D) 3 ]; (1) 

h/D = h 0 /D- 4(h 0 /D)\ 1.5- h 0 /£>)/3, 

h 0 /D=[l-(\-d 2 /D 2 f 5 ]/2; (2) 

S 0 =6(l-P){l-Z[h 0 /D- d/D(h 0 - h)/D]\/D- (3) 

q = 2/{P[\-(\ -2h/D)\ n P]/v)\; (4) 

r = 2P[(6/Zf - 2 /e 5 ]/S 0 ; (5) 

N c = 3Z (1 - P)/{ tlD 3 ),S oc = 3Z( 1 - P)d 2 /(4 D 3 ), (6) 

where h is the interparticle contact height; h/D is the normal 
deformation of the particles; v is the shape factor, equal to 
0.5 for spherical particles; h Q is the height of the spherical 
segment per contact. 

Calculations preformed using the expressions (1) and (2) 
with d= 0 agree with the experimental data. Thus, for 
Z=2.64 (percolation threshold) we obtain V= 0.157 while 
for Z = 7.3 we obtain the random packing density of the filler 
V r = 0.63. 

The structure of porous and composite monodispersed 
materials is often evaluated only according to the values of 
the average particle size of the disperse phase D , particle 
number density N (or particle volume density in a disperse 
medium (system) V = nD 3 N/6) and average distance L be- 
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tween the centers of the particles [6-10], which complicates 
the optimization of the structure of these systems. 

The simplest approach to determining the average inter¬ 
particle distance in a model of spherical particles is presented 
in [6]. If the relative volume V of the particles with diameter 
D in the system is known, then the following computational 
expression is obtained: 

D/L = (V/(n/6)f\ (7) 

If the number N of spherical particles per unit volume of 
the composite is known, then the calculation simplifies even 
more: L = 1/A 1 / 3 . 

It follows from the expression (7) that for V= V cr = 0.157 
we obtain D/L = 0.6693, for V=V r = 0.63 we have 
D/L= 1.0636, and only for V= 0.5235 the ratio will be 
D/L = 1.000 (V cr is the critical volume density of the parti¬ 
cles). In the experimental data, for example [5, 11], 
D/L = 1.000 in all cases for these values of V, i.e., the assess¬ 
ment obtained for the structure of the material is too arbi¬ 
trary. 

Attempts to refine the assessment of the structure of 
monodisperse systems in obtaining composites with pre¬ 
scribed properties and optimizing the structure of materials 
are examined in [7-10]. Thus, in [7-9] analysis of the 
structure of tribotechnical composite materials was used to 
develop durable powder materials with high antifriction 
properties. In [10] the problem of assessing the structure of 
monodisperse systems was analyzed from the standpoint of 
astrophysics. 

It is shown in [8] that the average distance between the 
centers of neighboring particles for a random uniform distri¬ 
bution of noninteracting particles, determined according to 
[10] as L = r(4/3)(47iA _ { / 3 ) = 0.5547V - ! / 3 , where N is the 
number density of particles, T(4/3) is the value of the 
gamma function for argument 4/3, as well as the ordinarily 
used estimate of the average distance according to the ex¬ 
pression (1) and a similar expression for hexagonal packing 
L = D = 1.2257V - / 3 are too rough. For this reason the fol¬ 
lowing interpolation relation was proposed: 

L(V) = a(V)N ] 0 n = [0.56 F 1/3 + 1/(2.24F 1/3 )] D, (8) 

where the function a(V) assumes the value 0.5554 as V —» 0 
and 1.1225 at F=0.74. 

The expression (8) was obtained by computer modeling. 
However, this algorithm is applicable only to density V= 0.45, 
showing its limited nature, since the actually achievable den¬ 
sity of random packing of noninteracting spherical particles 
V r is 0.63 [11]. 

In [10] the parameter L is calculated by a probabilistic 
method 

L= \.\9D/V l/ \ (9) 

where T(7/3) = 1/19 is the value of the gamma function for 
argument 7/3. 
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The use of the Poisson distribution makes the estimate of 
the parameter L less reliable. For example, for random pack¬ 
ing of a filler with V r = 0.63 the average distance between the 
centers of the spherical particles is equal to their diameter D. 
According to the expression (9) we obtain D/L = 0.7204. For 
critical filler density V cr = 0.157 we have D/L = 0.4534 and 
for V= 0.5 the ratio D/L = 0.6667. This also does not corre¬ 
spond to real systems [11]. 

The foregoing analysis shows that a more accurate 
method of evaluating the parameter L in composite systems 
must be developed. 

In the present work a more rigorous method of estimat¬ 
ing the parameter L in heterogeneous systems is presented. 
This method is based on the assumption that it is not suffi¬ 
cient to prescribe a single parameter L as the structural char¬ 
acteristic of a system. It is also necessary to give the coordi¬ 
nation number Z. 

Indeed, for volume density V of the disperse phase less 
than the random packing density V r the parameter L can de¬ 
pend on the volume density of the disperse phase and on the 
average coordination number Z (number of particles with 
which a central particle is in contact in the first coordination 
zone), and for composites also on fabrication technology, for 
example, on the thickness of the cover layer pre-deposited on 
a particle (if it consists of the matrix phase). Thus, in calcu¬ 
lating the relative volume of the solid phase of the system us¬ 
ing the expression (1) and (2) we obtain different values of V 
for a heterogeneous system for L = D : V=0.63 (Z= 7.3); 
V= 0.4 (Z = 4 ); V= 0.157 (Z= 2.64) with h/D = 0. 

For volume density V of the disperse phase greater than 
the random packing density V r the dependence of the param¬ 
eter L on the porosity is also not single-valued. The parame¬ 
ter L depends on the average size d of the interparticle con¬ 
tacts (or the normal deformation h/D at the center of the con¬ 
tact h/D), which for prescribed concentration of the solid 
phase depend on, in turn, the average coordination number Z. 
For definite types of vibratory compaction and d close to 
zero (the parameter L remains practically unchanged and 
equal to the average diameter of the particles D) the relative 
density of the system can be increased to V ~ 0.74 only as a 
result of an increase of the coordination number Z. 

From the standpoint of percolation theory low-concen¬ 
tration systems with density V aux < V cr = 0.157 (Z< 2.64) are 
unbound and hence indeterminate (F aux is the random pack¬ 
ing of auxiliary spheres with volume density of the particles 
with critical coordination number Z=2.64 and F aux = V CY = 
0.157). 

To increase the determinateness of the systems and to 
take account of the particle size and the type of system 
formed, we shall employ a system model in the form of 
quasipacking of auxiliary spheres with diameter L. It is con¬ 
venient to take as the base packing a random packing of 
auxiliary spheres with relative particle packing density V mx 
with critical coordination number Z=2.64 and F aux = V cr = 
0.157. Then the auxiliary spheres form a minimally bound 
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(statistically determinate) system along which conduction (of 
heat or electricity) is possible from the stand point of perco¬ 
lation theory. Since the average size D of the particles is less 
than Z, a partially determinate system of randomly distrib¬ 
uted noninteracting particles (each particle of system can be 
found anywhere inside the corresponding auxiliary sphere 
with diameter Z — on average at its center) is obtained. 
From the relation V aux = tiAZ 3 /6 = 0.157 with Z = 2.64 and 

V = tiAZ) 3 /6 we have 

L = 0.5395 DV~ 1/3 = 0.6693 N ~ 1/3 

forZ= 2.64 and F< 0.157. (10) 

The shortest average distance between the boundary of 
the particles 5 = L - D is given by the relation 

8 = (0.5395 F“ 1/3 - 1)£> = I[1 - (F/0.157) 1/3 ] 

forZ= 2.64 and F< 0.157. (11) 

The expressions (10) and (11) are also applicable for gas¬ 
eous or liquid finely disperse (colloidal) and course systems. 

For V aux = 0.157 with a cover layer of thickness 8 lay 
pre-deposited on the particles the critical volume density of 
the disperse phase V c ' r does not reach the theoretical critical 
value equal to 0.157 (see Fig. 1, dashed line): 

F c ' r =0.157(1-5 lay /Z) 3 . (12) 

It follows from the expressions (10)-(12) that for 

V < 0.157 the system is determinate only if the average size 
D of the particles or ( D + 5 lay ) equals Z and V aux = V cr = 
0.157. 

For F> 0.157 the system can be assumed to be statisti¬ 
cally determinate, i.e., consisting of randomly and uniformly 
distributed contacting particles. The value of Z is constant. It 
can be greater than D (if 8 lay > 1, i.e., the system being ana¬ 
lyzed is of the matrix type), equal to D (8 = 0, h/D = 0), or 
less than D (if h/D > 0). The coordination number is calcu¬ 
lated from the expression (l).ForZ=7.3 and V=V r the con¬ 
tacting particles form a rigidly bonded system and there is no 
longer any need to use a model of a system in the form of 
quasi-packing of auxiliary spheres with diameter Z. 

However, by analogy with the approach to loose packing 
it is convenient to examine the system being analyzed for 
V r > F> 0.157 as partially determinate, using the model of a 
system in the form of random close packing of auxiliary 
spheres with diameter Z with constant Z= 7.3. In this case 
only auxiliary spheres form a rigid bounded (determinate) 
system. The value of Z is also greater than D. The density of 
the system determines the shortest average distance 8 be¬ 
tween the boundaries of the particles: 

L = 0.8573W 1/3 = [3.78/(ttA0] 1/3 

forZ= 7.3 and 0.157 < F aux < V r ; (13) 

8 = (0.8573 V~ 1/3 - \)D = L[\ -(F/0.63) 1/3 ] 

for Z = 7.3 and 0.157 < F aux < V r . (14) 


The expressions (13) and (14) make the structural assess¬ 
ment of inhomogeneous systems more rigorous. 

For V=V r (in the presence or absence of a cover layer of 
thickness 8 lay pre-deposited on the particles) the particles of 
the solid phase form a rigidly bound (determinate) system. 
The value of Z can also be greater than D if the system being 
analyzed is of the matrix type (for example, matrix compos¬ 
ite). The shortest average distance 8 between the boundaries 
of the particles in such a system is determined by the thick¬ 
ness S lay of the cover layer pre-deposited on the particles. 

For F aux = V r and thickness 8 lay of the cover layer pre-de¬ 
posited on the particles the volume density of the random 
packing of the particles with a cover layer of a disperse phase 
V/ = V r (l - 8 lay /Z) 3 does not reach the value V r (see Fig. 1, 
dashed line). 

Subsequent packing (for example, vibratory compaction) 
of a disperse system, i.e., for V>V r , does not increase the 
statistical determinateness, with the exception of the extreme 
case (F= 1, Z= 14 [11]) where the expression (1) gives 
L/D = 1.0928 (see Fig. 1). Values of Z equal to or less than D 
and the porosity of the system are determined not only by the 
coordination number but also by the size of the interparticle 
contacts. Thus, the statistical determinateness of the system 
can be increased only if aside from Z and Z the average size 
of the interparticle contacts (or average normal deformation 
of the particles in the contact zone) is given. 

In the present model this can be taken into account by 
supplementing the expression (1) for porosity P < 0.37 by an 
approximate linear dependence on Z of the average normal 
deformation of the particles in the contact zone (we neglect 
the large difference in the sizes in the eight primary and six 
secondary contacts between the particles): 

h/D = (Z- 7.3)/158; 7.3 <Z< 0.14. (15) 

The maximum value Z= 14, adopted on the basis of the 
experimental data for compact material [11], substituted into 
the expressions (15) and (1) gives the relative density V= 1. 

Let us compare the value D/L = 1.0928 obtained from 
the expression (1) for V= 1 with the results following from 
the model of nonporous powder materials (Z = 14), based on 
elementary semi-regular polyhedral in the form of a trun¬ 
cated octahedron (bounded by eight regular hexagons and six 
squares). The polyhedral chosen completely fill the space. A 
coordination number of 14 and the average number of verti¬ 
ces in its faces, equal to 5.143, correspond to the experimen¬ 
tal data [11]. 

Nonporous material is obtained in this model for Z= 14 
with average deformation h/D = 0.04244 and correspond¬ 
ingly the contact size d = 0.5657D according to the expres¬ 
sion (1). It follows from the equality of the volumes of a 
semi-regular polyhedron and a sphere that for P = 0 we have 
h/D = 0.0395 or D/L= 1.0858. The values obtained are 
quite close, which supports our approach. 



Structural Assessment of Monodisperse Ceramic Composites 


115 



Fig. 1. Relative average distance between particles D/L in a disperse 
system versus the volume density V of the particles in this system. 

Similarly to the approach used to obtain the computa¬ 
tional expressions (13) and (14) for Z < 7.3 with system po¬ 
rosity less than the random packing the average distance L is 
best determined on the basis of the maximum coordination 
number Z = 14. 

To simplify the computational expressions we shall use 
the least-squares method for L with reference values of the 
relative particle packing density in the system V z=73 = 0.63, 
L Z=14 =L24D; V z= $ = 0.652, L z = 14 = 1.1 &D; U z=12 = 

0.7405, L z= 14 = 1.0905Z); V z= 14 = 1, L z= 14 = 0.9151. Since 
the average distance should not depend on the contact size 
(we are considering in each case a region where 14 particles 
are in contact with one particle), the contacts were assumed 
to point contacts (with the exception of the values for 
Z= 14). Thus, for V z = 12 = 0.7405 (closest hexagonal pack¬ 
ing) we have 12 direct point contacts (L = D) and two parti¬ 
cles, located from the chosen particle at distance 
L = 1.6337Z) (equal to twice the height of a regular tetrahe¬ 
dron). Then the average arithmetic distance is Z z = 12 = 
(12 + 2 x 1.6337)/14 = 1.0905. We obtain the result 

L = (2.25 - 1.335 0,5 ) D (for Z = 14), 

0.63 <V< 1 (L > 0.9151T>); (16) 

h = (0.6675 F° 5 - 0.625) D (forZ=14), 

0.63 <V< \ (h< 0.0424479). (17) 

Figure 1 shows the computational results for the average 
distance L between the centers of the particles according to 
the expressions (10), (13), and (16). The figure also indicates 
the experimental values V cr = 0.157 for Z= 2.64, V CY = 0.63 
for Z=7.3, and V=0.63 for Z= 12 (packing or vibratory 
compaction), obtained with different volume density of the 
filler for the same value of D/L = 1. In addition, the figure 
shows the results of the calculation of L from the expressions 
(7) and (9). The figure confirms the methods of [6 - 10], i.e., 
structural assessments of the monodisperse systems using 
only the average interparticle distance, are incorrect. 


CONCLUSIONS 

It has been shown that the conventional method used to 
assess the structure of monodisperse systems by means of the 
average interparticle distance is incorrect. This method was 
refined by indicating the average coordination number of a 
monodisperse powder system as well as by taking account of 
the connectedness of the system (V< V cr ) based on the cor¬ 
puscular model and percolation theory. 

A refined method of determining the average distance 
between particles in dense disperse systems for V > V r taking 
account of the coordination number and the average size of 
the interparticle contacts was proposed. 
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